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ABSTRACT 

With improving methods and surveys, the young field of extrasolar planets has 
recently expanded into a qualitatively new domain - terrestrial (mostly rocky) 
planets. The first such planets were discovered during the past year, judging 
by their measured masses of less than 10 Earth-masses (M®) or Super-Earths. 
They are introducing a novel physical regime that has not been explored before as 
such planets do not exist in our Solar System. Their composition can be either 
completely terrestrial or harbour an extensive ocean (water and ices) above a 
rocky core. We model the structure and properties of the first Super-Earth 
(mass ~ 7.5 M®) discovered in 2005, illustrating the possibilities in composition 
and providing radius evaluations in view of future detection of similar planets 
by transits. We find that a threshold in radius exists for which larger values 
indicate that a Super-Earth most certainly has an extensive water content. In 
the case of GJ876d this threshold is at about 12000 km. Our results show that 
unique characterization of the bulk composition of Super-Earths will be possible 
in future transit studies. 

Subject headings: planetary systems — planets and satellites: individual GJ876d 
— Earth 
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1. Introduction 

In the last 11 years more than 200 exoplanets have been discovered but only in the last 
year the minimum mass for detection has been pushed down to allow for terrestrial planets. 
The first extra-solar planet with a mass lower than 10 M e , i.e. a Super-Earth, was discovered 
last year by Rivera et al. (2005) orbiting an M4V star named GJ876. Its estimated mass 
is 7.5±0.7 M e and it has an orbital period of 1.94 days. It is close to the host star and 
the surface temperature is calculated to range between 430 and 650 K (Rivera et al. 2005) 
. The star is known to have two gas giant planets GJ876b,c in 30 and 61 day orbits (Marcy 
et al. 1998, 2001). Hence, the planetary system of star GJ876 has both an architecure and 
a planet - the Super-Earth, that are unfamiliar to our Solar System experience. Two other 
low mass planets have been discovered in this year: OGLE-2005-BLG-390Lb (~ 5M® at 5 
AU — (Beaulieu et al. 2006)) and HD69830b (~ 10M e at 0.08 AU — Lovis et al. (2006) 
) exemplifing the variety in Super-Earths that will be discovered in the near future with 
missions like Kepler. 

Calculating the internal structure of a Super-Earth can help us determine how different 
or similar to the Earth this planet might be. One of our goals is to provide mass-radius re- 
lations that will help characterize Super-Earths discovered by the transit method. We show 
in this study the tools we have to characterize Super-Earths and use GJ876d as an exam- 
ple of our model capabilities and limitations. We first describe the numerical method used 
to obtain density, pressure, and temperature profiles as a function of radius and refer the 
reader for more details to Valencia et al. (2006), where we provided the first theoretical grid 
of models for Super-Earth planets. To calculate the planet's radius and internal structure we 
need to make reasonable assumptions about its composition. It might be completely rocky 
or might have accumulated a substantial amount of ices depending on the material available 
during formation. We show the results and implications for different likely compositions for 
GJ876b as an example that can be replicated for any Super-Earth. Additionally, we explore 
the effects of tidal heating that might be present in short period orbits (like GJ876d) and 
nonzero eccentricities. 



2. Model 

2.1. Numerical method and Equations of State 

We model the planet as being composed of distinct, homogeneous in composition, spher- 
ical shells. For a terrestrial planet these shells are: the mantle that is divided into lower 
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mantle and upper mantle, and the core that depending on the temperature structure of the 
planet might be divided into a liquid outer core and a solid inner core, as is the case for the 
Earth. The composition for the different layers is taken from the mineralogical composition 
known for the Earth (McDonough & Sun 1995). The upper mantle is composed of olivine 
(ol) and higher pressure forms of olivine (wadsleyite (wd) and ringwoodite (rw)); the lower 
mantle develops when rw transforms to perovskite (pv) and ferromagnesiowustite (fmw), 
with an additional shell at high pressures when pv transforms to post-perovskite (ppv). The 
solid inner core in the Earth is composed of Fe and small quantities of Ni, and the outer 
liquid core is composed of Fe and a light alloy (see figure 1). The candidates for this alloy are 
S, Si, O, C and H and to this day there is no consensus on which one(s) or the amount(s). We 
have used in this study a composition of pure Fe and Feo.s(FeS)o.2 to show the uncertainties 
in the radius of a planet from the lack of knowledge of the composition of the core. 



For planets that harbour a substantial H20/ice content (ocean planets), there is an 
additional layer overlying the rocky interior (mantle and core), composed of two shells - 
water above high pressure phases of ice. In some cases pressure is not high enough (amount 
of H 2 is not large enough) to solidify water and the ice layer is absent (figure 1). The 
thickness of the water is determined by the intersection between the PT curve of the planet 
and the melting curve of ice. In the case of low surface temperature, the planet will have 
an additional shell above the water layer composed of ice I — the lightest form of H 2 and 
negative Clapeyron slope - as is the case for Jupiter's satellites. 



The numerical model solves the following differential equations for density p, gravity g, 
mass m and pressure P in each shell with a Runge-Kutta 4th order solver: 



dp p{r)g{r) ^ 



dr 0(r) 
dg A „ . . 2Gm(r) 
dr r 



A*Gp(r) - (2) 



^ = 4irr 2 p(r) (3) 
or 

dP 

— = -p(r)g(r) (4) 



where 0(r) = is the seismic parameter that can be calculated with an equation of 
state (EOS) for K S: the adiabatic bulk modulus; G is the gravitational constant and r is the 
radius. 
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The model integrates from the surface inwards with boundary conditions for P, T, p 
and total mass (M). It starts with a guess for the planet's radius (R) that yields a surface 
gravity value and determines the structure to the last shell where there could be excess or 
deficit in mass depending whether the starting R is too small or too big respectively. The 
model then uses a bisect-newton approach to determine the radius that yields zero mass and 
zero gravity at the center of the planet (within 5 meters). 

The EOS that we have implemented for all layers except the water layer is the Vinet 
EOS, because is the best analytical EOS for extrapolation (Hama & Suito 1996): 



K T (p, 300) = 3K {x 2/3 - x 1 / 3 ) exp - 1) (l - x' 1 ^ (5) 

where x — — , Kn and K' n are the isothermal bulk modulus and its first derivative at a 
reference state — zero pressure and 300 Kelvin — and K s = K T (1 + eryT), where a is the 
thermal expansion coefficient and the Gruneisen parameter 7 = 7o:r~ 9 with q = — Most 
high-pressure experiments fit their data to the 3rd order Birch-Murnaghan EOS (BM3) that 
is derived from a truncation of the linear expansion of the Helmholtz free energy with strain 
to the third order (Poirier 2000). The extrapolation of this EOS is highly uncertain since 
the fourth term in the expansion might not be smaller than the 3rd term in some cases. We 
use the values for K and K obtained from the literature corresponding to the BM3 EOS 
to obtain P(V) and refit the data to the Vinet EOS within a volume compression between 1 
and 2/3. The general trend is that K Q is smaller by a few GPa and K' Q is larger by ~6-8% in 
the Vinet fit. We used the Rangkine-Hugoniot EOS from Stewart & Ahrens (2005) within 
the water region. 



To incorporate the effects of temperature in the EOS we added a thermal pressure 
term (due solely to deviations in temperature from 300K) for core and mantle regions that 
translates into a thermal bulk modulus correction. With R the universal gas constant, 
9 = 9q exp ( \ as the Debye temperature and n as the number of atoms in the unit cell, 



the thermal correction is expressed as: 



K T (p,T) = K T (p,300) + AK th (p,T) (6) 

AK th = 3nR 1P (f(T)-f(T )) (7) 

/(T) = (l-g-3 7 )- / ^-^- + 3fl 7 (8) 

v 6 Jo exp £ — 1 exp (v/T) — 1 



For the icy region we incorporated the temperature effects in the density via the thermal 
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expansion coefficient: 



p(P,T) = p(P,300)exp / a(P,T')dT 



(9) 



J 300 



a(P,T) = (ao + aiT) 1 + -^P 




(10) 



where a , ai and 6 are coefficients determined at zero pressure. 



2.2. 



Thermal model 



In order to incorporate the temperature effect in the EOS we need to have a model 
describing the temperature profile in the planet. We use the Earth as our starting point and 
parameterized convection theory to model the temperature regime of Super- Earths. The sur- 
face heat flow on Earth (44TW Pollack et al. (1993))reflects the heat from radioactive sources 
in the mantle (some authors also place potassium in the core) and secular cooling from an 
initial hot state. If we account for the concentration of heat sources in the crust (generating 
about 5-12TW (Davies 2001; Vacquier 1991; Taylor & McClennan 1998), we take 8TW), and 
assume a bulk silicate concentration of uranium, potassium and thorium, radioactive heat in 
the mantle accounts for ~ 58% of the total heat flow. In Valencia et al. (2006) we reported 
scaling the heat flow with mass for Earth-like planets as a first attempt to scale radioactive 
heat sources and secular cooling together. Here we refine this scaling by separating the 
treatment of the two sources. If we assume that the material that makes Super-Earths is the 
same as the bulk material for the Earth, then the radioactive heat sources scale with mantle 
mass. Planets that have larger mantles will have higher heat flows. It is difficult to scale 
secular cooling for a Super-Earth because we do not know the role that mass plays in the evo- 
lution of a planet. With the scaling laws derived in Valencia et al. (2006) and parameterized 
convection, we have found (Valencia et al, in prep) that a massive planet is likely to convect 
in a plate tectonic regime similar to the Earth. Intuitively, the more massive the planet 
is, the higher the Rayleigh number that controls convection, the thinner the top boundary 
layer (lithosphere) and faster the convective velocities. This scenario might aid subduction 
of the lithosphere causing the onset of plate tectonics. Thus, we adopt that the evolution 
of a Super-Earth might be one that leads to the same proportion of secular cooling to ra- 
dioactive heating as in the Earth. This is the first attempt to scale secular cooling with mass. 



We define the Rayleigh number Ra (dimensionless parameter that controls convection) 
in terms of the heat flux (q s ): 
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Ra = 



pgaq s /k 



(11) 



where g, k, k, 77 are the average gravity, conductivity, diffusivity and viscosity of the 
mantle. The thickness of the top boundary layer is 



independent of the size of the mantle, where Ra crit ~ 1000 is the critical Rayleigh number 
and a is a coefficent of order unity. We treat viscosity in two ways: i) isoviscous case and 
ii) a temperature dependent treatment r](T) = r] (jr)~ 30 (see Valencia et al. (2006) for more 
details). 

The planet modeled in this study has a convective core, mantle and ice/water layer 
where the temperature can be described as adiabatic: 



All interfaces between the chemically distinct layers develop boundary layers that trans- 
fer heat conductively 

%-\ <*> 

Their thickness depends strongly on the local viscosity. Owing to the low viscosity of 
the core and water we only model conductive boundary layers at the top and bottom of the 
mantle. The top boundary layer thickness is described in equation (12). In an isoviscous 
internally heated case, the bottom boundary layer of the mantle would vanish. Conversely in 
an isoviscous heated-from below case it would have the same thickness as the top one. Since 
we assume these planets are in a intermediate regime (same as Earth) the bottom boundary 
layer is taken to be half of the top one, and the heat flux reflects the heat coming out of the 
core. 

Even though there are two critical assumptions in determining the temperature structure 
- the constant ratio between internal heating and secular cooling, and the thickness of the 
lower boundary layer in the mantle — the effects of temperature in the internal structure 
of a planet are small, particularily the effects in radius. The temperature profile is needed 
mostly to determine the location of the phase changes that determine the different regions 
in the planet. A temperature dependent viscosity increases the radius by ~ 20 km compared 
to the iso- viscous case showing that neither assumption is critical. 




(12) 




(13) 



-7- 



The numerical model iterates several times until convergence exists between the bound- 
ary layer thickness and the radius (i.e. surface heat flux) so that the temperature profile is 
determined self-consistently. 

The numerical model needs as input the composition of the different regions as well as 
the proportion of ice mass fraction (IMF) and core mass fraction (CMF), and surface P and 
T. The output is: density, pressure, temperature, mass and gravity as a function of radius; 
the total radius and the location of the phase changes and composition boundaries. Table la 
shows the different values in composition used in this study. The composition assumes that 
the mantle minerals have 10% of iron and 90% magnesium in the mineral structures, except 
for post-perovskite phase that is believed to incorporate more Fe (Mao et al. 2004) where we 
adopt a 20% Fe by mol composition. Table 2 shows the values we used to determine the phase 
boundaries between the major silicate phase changes (ol to wd, rd to pv+fmw, pv to ppv), 
the phase boundary between water and ice is determined from a Krauss-Kennedy melting 
curve proposed by Stewart & Ahrens (2005). We calculate the melting point of iron by using 
a Lindeman equation (see Valencia et al. (2006) for details) only to test whether the core is 
in a solid or liquid phase but not to distinguish between different phases within the core (as 
is the case for the Earth). Only the Earth-like composition yields a planet hot enough to 
have a liquid core (see figure 3). A planet with more water content (IMF=40%) will have 
a smaller mantle for a given CMF and will be cooler (due to less radioactive heating) than 
a less water rich planet (IMF=20%). The strong dependency of the Gruneisen parameter 
with volume in the mantle (pv or ppv) (that determines the adiabatic gradient) prevents the 
temperature at the core-mantle boundary of being hot enough in planets with medium to 
small mantles (ocean planets with medium cores and rocky planets with large cores). The 
adiabatic gradient within the mantle remains small and the temperature increase within the 
mantle is not sufficient to allow for a liquid core (figure 3). 



2.3. Composition cases for GJ876d 

We consider a few basic scenarios for the composition of GJ876d that are representative 
but probably not exhaustive. They are (1) a simple Earth-like composition, a planet with a 
core that makes up 33% of the total mass, with a mantle that is 10% (by mol) iron enriched, 
and a lower mantle composed of 30% fmw and the rest pv. (2) an Earth-like, but with a 
very large core (mass fraction of 80%), (3) an ocean planet with a 20% by mass water/ice 
layer on top of a terrestrial (Earth-like) core, and (4) an ocean planet with a 40% water/ice 
layer. Our reasons for introducing the later three scenarios are: case (2) is analogous to 
Mercury - with its very close orbit, GJ876d might have similarly acquired a massive core 



- 8- 



during its formation - this case is essentially a Super-Mercury (see Valencia et al. (2006) for 
details). Cases (3) and (4) are inspired by the possibility (Zhou et al. 2005) that GJ876d 
actually formed beyond the snow line in the protoplanetary disk of GJ876. The snow line 
in a disk is the distance from the star at which the local temperature and pressure in the 
midplane allow ices to exist during planet formation (Sasselov & Lecar 2000). In our Solar 
System the snow line was at about 2-3 AU, but GJ876 is a small 0.32 MO star (Rivera et al. 
2005). With its lower energy output and smaller size, GJ876 would define a snow line at 
~0.25 AU in its passive disk at zero age. Beyond the snow line, ices (mostly H 2 0) increase 
the solid fraction in the disk by a factor of ~ 4, given the solar composition of the star. 
Therefore a Super-Earth forming in this region of the disk could acquire and retain water 
as a major fraction (20% or more) of its total mass. Differentiation will quickly produce an 
Earth-like core and interior overlaid by water ice (higher pressure ice phases such as VII and 
X) and a deep water ocean, similar to the ocean planets described by Leger et al. (2004). A 
more complicated issue is whether this water could evaporate during the lifetime (~6 Gyr) 
of the planet. A simple thermal escape estimate precludes that. It is possible to retain an 
icy/liquid layer over the age of this planetary system with a surface temperature of ~550 
K provided the atmospheric pressure is large enough. Ice VII is stable at such temperature 
given pressures larger than a few GPa (Stewart & Ahrens 2005) . Furthermore, the atmo- 
spheric partial pressure would have to exceed ~100 MPa to prevent the water layer from 
evaporating (Wagner & Pruss 2002). A non-thermal escape calculation is beyond the scope 
of this paper; it could prove very efficient in removing a lot of water from GJ876d, especially 
during the early evolution of the star. 

Addressing the issue of bulk planet composition is central to our study. The case of 
GJ876 provides us with useful constraints. First, we know (Rivera et al. 2005) that the star 
has a composition very similar to the Sun and is of similar age. Second, there is evidence that 
the orbits of the GJ876 planets evolved shortly after planet formation by inward migration 
- the very close orbits of all 3 planets and the 1:2 mean- motion resonance state of the giant 
planets (Zhou et al. 2005) 



3. Results 

For the Earth-like composition case (lines with stars in top figure 2) perovskite trans- 
forms to post-perovskite at a radius of ~ 9800 km, so most of its mantle is actually composed 
of ppv+fmw. This brings attention to this relatively newly discovered phase present in the 
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lower-most mantle on Earth (only between 125 to 136 GPa) and the importance of determin- 
ing with accuracy its behaviour at high pressures. If this planet had a pure Fe core, its radius 
would extend to 10800 km or 1.70 i?© (Earth radius) and an additional 130 km (1.71.R©) 
if the core was composed of Feo.s(FeS)o.2- Despite the large differences in core density (300 
kg/m 3 ) between the two core compositions, due to the different locations of the core- mantle 
boundary that satisfy the CMF, the radius differs by very little. This implies that the radius 
for planets with medium sized cores is a robust parameter. 

The second scenario looks at the effect of having a larger CMF of 80% (see top figure 
1) where we expect the composition of the core to have a larger effect. A planet that has a 
larger iron content presumable mostly in its core, can accomodate more of its mass in the 
central region and therefore have a smaller radius compared to an iron-poor planet. This 
type of planets (Super- Mercuries, Valencia et al. (2006)) exhibit comparatively large bulk 
densities as Mercury does in our Solar System. If GJ876d had a core of 80% by mass its 
radius would only be 9200 km (or 1.45 R 9 ) for pure Fe core and 9600 km (or 1.5 i?©) for a 
sulfur enriched core. 

Next consider the possibility of this planet having a water/ice layer amounting to 20% 
and 40% of the total mass on top of a terrestrial core (bottom figure 2). Despite ice and water 
being highly compressible (bulk modulus is 1-2 orders of magnitude smaller respectively than 
for silicates), their densities at initial compression are low and the relation of the density 
gradient to density is quadratic. Consequently, a large amount of water makes the planet 
larger. If the Earth had a substantial amount of water the radius would extend to 7100 
km (1.11-R©) for 20% water content and 7600 km (1.19/?©) for 40% water content. Figure 2 
shows the density structure and figure 3 - the pressure-temperature structure for the different 
composition cases considered above. The ocean planets have a layer of water above a layer 
of denser ice due to the positive Clapeyron slope of ice VII (note the melting curve of ice VII 
in figure 3). Hemley et al. (1987) showed that ice transforms gradually from ice VII to ice X 
with pressure and proposed an EOS for the combined system. We have adopted those values 
for this calculation. Below the ice layer the pressure is ~ 150 GPa, beyond the transition 
pressure between pv and ppv (P=125 GPa at T=2750K Tsuchiya et al. (2004)), therefore 
the silicate mantle is made of ppv+fmw — perovskite and upper mantle silicate phases are 
absent. If GJ876d had a water/ice layer that accounted for 20% and a core that accounted 
for 33% of the total mass, the radius would be 11900 km (or 1.87 R e ). An additional 20% of 
water would extend the radius 550 km more (to 1.95 it!©). These radii are 11 and 15% larger 
respectively, compared to the Earth-like composition radius result. Such planets will be eas- 
ier to detect when transiting their stars. Adding sulfur to the core would extend the radius of 
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both ocean planets by ~130 km. The radius results for all compositions are shown in Table 3. 

In summary, the four bulk compositions produce differences in the Super-Earth mean 
densities and radii that are large enough - from 3% to 30% in radius, to be observationally 
measurable. For example, if GJ876d was transiting its star, MOST satellite photometry 
(Rowe et al. 2006) could achieve 3% precision in its radius, with uncertainty dominated by 
the stellar radius (R = 0.32R Q ). 

Our main result is that with a measurement of both mass and radius (even with a pre- 
cision of 10%) we could certainly tell between an iron rich differentiated terrestrial planet 
(R < 9600 km) and a water rich (IMF>20%) ocean planet (R > 12000 km). This is because 
there is a limited family of compositions that would fit both. To illustrate this point we 
have calculated the trade-off curve between the amount of ice (IMF) and core (CMF) that 
would satisfy M=7.5M e and R=11900km (see Figure 3). The maximum of amount ices this 
planet could harbour is ~ 49% with CMF=51% and no mantle . The minimum could be 
zero if the planet had no core. This extreme case (no core) is so unlikely that a radius detec- 
tion of 11900km or more would necessarily mean that the planet has a water layer. Hence, 
knowledge on just the radius and mass of a planet - in combination with internal structure 
models — could yield valuable structural and compositional information on the planet and 
its formation environment. 



The Super-Earth's extreme close proximity to its star GJ876 means that even small 
perturbations on its orbit might induce significant additional heating from tides and change 
its internal structure. We check this by a simple analysis, with a surface heat flux consisting 
of energy from radioactive heat sources, secular cooling and tidal heating. We add different 
amounts of tidal heating Qudai to the heat flow Q (see section 2.3) that translate in higher 
Rayleigh numbers and thinner lithosperes and calculate the thermal structure: 



For the calculation of the effects of internal heating we assume a temperature dependent 
viscosity. A simple analysis states that the energy dissipation rate produced by an eccentric 
orbit in a synchronous planet is (Murray & Dermott 1999) 



4. Tidal heating effects 



Qtotal — Q + Qtidal 



(15) 
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where G is the gravitational constant, is the mass of the star, a is the semimajor axis, 
n is the mean motion, jx is the ratio of elastic to gravitational forces, jl ~ (10 4 km/i? p ) 2 , and 
Qd is the specific dissipation parameter. In other words the energy available for tidal heating 
of the 7.5 M e Earth-composition planet GJ876d is 

1.2 x 10 2 V [W] (16) 

where Qd = 280 as for the Earth (Ray et al. 2001). Figure 5 shows the effects of differ- 
ent values of dE/dt on the PT regime for this planet. In addition, figure 5 also shows the 
solidus for MgO, the lower mantle mineral with the lowest melting temperature. By in- 
creasing the amount of tidal heating the internal temperature rises and in some cases the 
lowermost mantle might be partially molten. Our calculations show that large increments 
in tidal heating cause small increments in internal temperature. Thus, in order to melt the 
lowermost mantle, the tidal heating needs to exceed ~ 6.8 x 10 17 W for this 7.5 M e planet (or 
~ 10 17 W/M e ), meaning 2500 times the heat flow without tidal heating. Even though the 
internal temperature is much higher than without tidal heating, the radius only increases by 
~100 km. For planet GJ876d melting could commence for eccentricites above 0.008. While 
we do not yet know the dynamic environment of this planet very well, such eccentricities are 
not excluded (Rivera et al. 2005). We note that the newly discovered ~10 M e Super-Earth 
HD69830b (Lovis et al. 2006) orbiting at 0.08 AU from a 0.86M© star is also in a 3-planet 
system and appears to have a non-zero e = 0.10 ±0.04. Hence its tidal heating could exceed 
10 17 W/M @ . HD69830b is expected to be mainly rocky (Lovis et al. 2006) and we estimate 
its radius at 1.84 i? ffi (Valencia et al. 2006). 



5. Uncertainties 

The uncertainties in the model come from uncertanties in the EOS and the composition 
of the planet that needs to be known at high pressures and temperatures. Any phase change 
has to be known a priori in order to be incorporated in the model. Therefore when we 
extrapolate to high pressures it is possible that the material might change into a different 
phase that we can not account for without experimental evidence. If there are high-pressure 
phases unaccounted for in the model, the radius here obtained for different compositions 
would be an overestimate owing to the denser character of high-pressure minerals. Benoit 
et al. (1996) describes a transition of H 2 to ice XI at 300 GPa. GJ876d with IMF=40% 
reaches a maximum pressure of 350 GPa at the bottom of the ice shell. Therefore, the effects 
of ice XI would be to decrease the radius by a few tens of kilometers for this planet and more 
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for planets with larger masses and IMF. 

The likely compositions chosen for GJ876d yield a central pressure between 2500 and 
5000 GPa. Values that exceed any laboratory experiment up to date (Cohen et al. 2000). 
The extrapolation of EOS — Vinet and 3rd order Birch-Murhnaghan — is yet to be tested 
at these pressures. Nevertheless, the uncertainty in composition is greater than errors in 
this extrapolation. The difference in radius from using the two different EOS is ~ 100 
km, less than the difference due to the different compositions considered here. There are 
also uncertainties in the determination of the parameters in the EOS that translate into 
uncertainties in the internal structure. We considered two different data sets for the EOS of 
pure Fe (Uchida et al. 2001; Mao et al. 1990) and find that the difference in radius is only 
of ~65 km for the planet dominanted by core composition. 

Finally, the different treatments in viscosity (isoviscous vs. temperature dependent) only 
yield a difference in radius of ~20 km, though the temperature regime is different for the 
four compositions described herein. Our lack of knowledge on the temperature structure is 
offset by the small temperature effects in the total radius of the planet. This means that 
despite the uncertainties in the thermal model (viscosity assumptions and heat flow scaling 
assumptions) the radius is a robust parameter. Therefore we place an estimate on the error 
in radius of ~200 km due to the lack of knowledge associated with the EOS and in the 
thermal profile. 

6. Conclusions and Discussion 

In conclusion, the new surveys and improved detection techniques have opened up the 
study of a new class of objects, Super-Earths, that brings us closer to characterizing the 
exoplanets most similar to Earth. In particular, we can model the internal structure of the 
first discovered Super-Earth by looking at likely compositions. We find that the radius varies 
between ~9200 and 12500 km. If the orbital geometry allowed a transit follow up to GJ876d, 
the expected flux drop would be of 1.9 — 3.6 x 10~ 3 , large enough to be observed. If observed 
with a large-aperture telescope (Holman et al. 2005) or from space, e.g. with MOST (Rowe 
et al. 2006), a transiting planet like GJ876d would yield a radius determination with ~ 3% 
precision and one would be able to distinguish between all four scenarios/models presented 
by us here, especially between an iron rich differentiated and an ocean planet. Moreover, 
for a given mass there is a radius that delimits the boundary between ocean planets and 
terrestrial planets. If GJ876d radius was larger than 12000km, it would indicate it is an 
ocean planet. 
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With the upcoming space mission Kepler that ability will improve, especially with ex- 
pected advances in RV measurements (hence planet masses) for Kepler's fainter targets. 
Therefore, when radii measurements of terrestrial planets become available, we can put con- 
straints on the bulk composition of a planet (given its mass) and begin to understand the 
conditions on distant terrestrial planets. A following article will treat this problem by show- 
ing in a simple manner the relation between radius, mass and composition. 
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Table 1. Compositional data 



layer 


composition 


po (kg/m 3 ) 


K (GPa) 


K 


7o 


Q 


#0 


ref 


OCeaii 




one oq 
yyo.zo 


O 1Q 
Z.lO 










i 
i 


shell 


ice VII+X 


1463 


2.308 


4.532 


1.2 


I. 


1470 3 


2,3 












ai=-4.2xl0~4 


a 2 =1.56 xlO-6 


b = 1.1 


2 


upper 


ol 


3347. 


126.8 


4.274 


0.99 


2.1 


809 


4 


mantle 


wd+rw 


3644. 


174.5 


4.274 


1.20 


2.0 


908b 


4 


lower 


pv+fmw 


4152. 


223.6 


4.274 


1.48 


1.4 


1070 


4 


mantle 


ppv+fmw 


4270. 


233.6 


4.524 


1.68 


2.2 


1100 


5 


core 


Fe 


8300 


160.2 


5.82 


1.36 


0.91 


998 


6,7 




Feo. 8 (FeS) .2 


7171 


150.2 


5.675 


1.36 


0.91 


998 


6,7 



References. — Compositional data used in each shell of the model has been taken from different sources 
and refit to the Vinet EOS: (1) Stewart & Ahrens (2005) — Rankine-Hugoniot EOS (2) Hemley et al. 
(1987) (3) Fei et al. (1993) (4) Stixrude & Lithgow-Bertollini (2005) — A reuss average was performed 
according to a mixture of 10% Fe and 30% fmw for the bulk modulus. The thermodynamic parameters 
(7O, q and 6q) were taken from the most dominant phase in each shell. (5) Tsuchiya et al. (2004) - an 
increase in density was needed to account for 20% Fe in ppv according to dln/dx = 0.3 where x is the iron 
content (Mao et al. 2004). A Reuss average was then used with fmw. (6) Williams & Knittle (1997) — we 
also tried Uchida et al. (2001) parameters for density and bulk modulus and found a discrepancy in the 
radius of only 60 km for a compositon of CMF=80%. (7) Uchida et al. (2001) 
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Table 2. Phase boundaries in silicate mantle 



transition 


boundary 




ol — > wd+rw 


T=400P-4287 




rw — > pv+fmw 


P=22.6 


if T > 1750K 




T=13573-500P 


if T < 1750K 


pv+fmw — > ppv+fmw 


T=133P-1392 





Note. - Table showing the different phase boundaries 
between the silicate mantle minerals used in this study. P is 
in GPa and T is in K. 



Table 3. GJ875d radius 



RADIUS 


Fe 


Fe . 8 (FeS)o. 2 


Earth- like 


10786 km 


10914 km 


Rocky CMF=80% 


9228 km 


9580 km 


Ocean IMF=20% 


11890 km 


12014 km 


Ocean IMF=40% 


12448 km 


12576 km 



Note. — Radius results for the different four 
bulk compositions for a pure Fe core and a core 
with Fe and S. The only appreciable difference 
of 400 km happens in a core-dominated planet. 
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Upper Mantle: olivine 

Lower Mantle: perovskite 
+ ferromagnesiowustite 



Outer core: 
Fe or Fe . 8 (FeS)o.2 

Inner core: 
Fe orFe . 8 (FeS) .2 



Lower Mantle: post-perovskite 
+ ferromagnesiowustite 

Upper mantle: 
wadsleyite+ringwoodite 




Fig. 1. — Schematic representation of the model. To calculate the internal structure of Super-Earths we 
assume a similar composition to the Earth's: A dense core of pure Fe or Feo.s(FeS)o.2 as end member cases 
(the Earth has an outer core of Fe plus some unknown alloy and the solid inner core has Fe and Ni); a lower 
mantle composed of two silicate shells (ppv-hfmw, pv+fmw); an upper mantle composed of two silicate shells 
(wd+rw, ol). The thickness of the shells will depend on the PT profile for the planet and the amount of 
mass in the core. An ocean planet - right - will have an additional water/ice layer above the rocky core. 
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Fig. 2. — Internal Structure of exoplanet Gliese 876d: Density Profile. Four different compositions are 
illustrated. The surface is on the right and the centre of the planet is on the left. The solid lines show 
the cases for which the composition of the core is taken to be pure Fe. Dashed lines are for the case of 
Feo.2(FeS)o.8- Lines with stars shows the internal structure of GJ876d if the composition was Earth-like. 
Square symbols show the density profile if this planet had 80% of the mass in the core. Circles show the 
structure if this planet had formed outside the snow line and retained 20% of its mass as a water/ice layer. 
Diamonds shows the density structure if GJ876d had retained 40% of water/ice. A Preliminary Reference 
model for Earth (PREM - Dziewonski & Anderson (1981) ) is shown for reference. 
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Fig. 3. — Internal Structure of exoplanet Gliese 876d: Pressure- Temperature Structure. The PT structure 
of GJ876d is illustrated considering the four compositions presented in this study and pure Fe in the core. 
The ocean-compositions cross the ice VII solidus at ^55 GPa which leaves means the water layer is 1200 km 
deep for both ocean planets (IMF=20% and IMF=40%). 
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Fig. 4. — Trade-off curve beetween CMF and IMF. For a planet with M = 7.5M® and R = 11900 km the 
trade-off curve between the amount of ices and amount of mass in the core is depicted as a solid line. The 
squares correspond to the different data points. The dashed line shows the amount of mass in the mantle. 
The maximum amount of ices this planet could harbour is 48.6% with CMF=51.4% and no mantle. A planet 
this size can only be rocky if it had no core (IMF=0 when CMF=0). 
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Fig. 5. — Effects of tidal heating on PT structure. PT structure for a 7.5 M e planet experiencing different 
degrees of tidal heating. Dots: no tidal heating. Squares: tidal heating that is 10 times the heat flow 
withouth tidal heating (Q) or total heat flow Q total =3000 TW. Downward-pointing triangles: 200 times 
or Q totaJ =5500 TW. Upward- pointing triangles: 2500 times or Q tota; =68000 TW. The solidus for MgO is 
crossed in the mantle when the tidal heating is 2500 times the radioactive and secular cooling contributions. 
Dotted lines refer to a thermal model without a lower boundary layer in the mantle. As the tidal heating 
increases the lower boundary layer becomes negligible. 



